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Graphene has attracted much interest in both academia and industry. The challenge of making it 
semiconducting is crucial for applications in electronic devices. A promising approach is to reduce its 
physical size down to the nanometer scale. Here, we present the surface- assisted bottom-up fabrication of 
atomically precise armchair graphene nanoribbons (AGNRs) with predefined widths, namely 7-, 14- and 
21-AGNRs, on Ag(lll) as well as their spatially resolved width-dependent electronic structures. STM/STS 
measurements reveal their associated electron scattering patterns and the energy gaps over 1 eV. The 
mechanism to form such AGNRs is addressed based on the observed intermediate products. Our results 
provide new insights into the local properties of AGNRs, and have implications for the understanding of 
their electrical properties and potential applications. 

Graphene, a single sheet of sp 2 -bonded carbon atoms arranged in a honeycomb lattice, exhibits extraord- 
inary electrical and mechanical properties. For applications in electronic devices, graphene needs to be 
semiconducting. Reducing its physical size down to the nanometer scale is a promising approach. 
Theoretically, free-standing armchair graphene nanoribbons (AGNRs) can behave as either metals or semi- 
conductors, depending on their widths 1 " 3 . GNRs fabricated by top-down methods generally possess uncontrolled 
and atomically imperfect edges, precluding detailed studies on finite size and edge effects 4-6 . Recently, the 
fabrication of 7-AGNRs (Fig. lb, the number before "AGNR" indicates its width, that is, the number of car- 
bon-pair lines in the direction perpendicular to AGNR axis, as sampled in Fig. Id) from 10,10'-dibromo-9,9'- 
bianthryl (DBBA, Fig. la) via polyanthracene has been demonstrated 7,8 . The GNR topology, width and edge 
periphery are predefined by the DBBA precursor molecule. Density functional theory (DFT) calculations suggest 
that dehydrogenation takes place in a domino -like fashion and is more favourable at sites between anthracene 
monomers than at their head sites 9 . This limits the formation of wider AGNRs such as 14- and 21-AGNRs (Fig. lc 
and Id), which can be regarded as two and three 7-AGNRs conjugated side-by-side via dehydrogenation at the 
head sites, respectively. 

Fabrication of atomically precise AGNRs with varying widths provides a unique opportunity to investigate 
quantum confinement and edge effects on them. In this article, we present the surface-assisted bottom-up 
fabrication of a series of atomically precise AGNRs from DBBA on Ag(lll), such as 7-, 14- and 21-AGNRs, 
as well as their spatially resolved width- dependent electronic structures. Their typical vibration modes are 
measured by Raman spectroscopy. The associated electron scattering patterns (ESPs) which arise from intervalley 
electron backscattering at both armchair edges 10,11 are visualized by scanning tunneling microscopy (STM). The 
local density of states (LDOS) on each AGNRs observed by scanning tunneling spectroscopy (STS) reveal the 
presence of energy gap over 1 eV. The mechanism to form such AGNRs is addressed based on the observed 
intermediate products. Our results are a significant step towards understanding the local properties of AGNRs, 
and provide us with design guidelines for the fabrication of edge-precise graphene nanoribbons for potential 
device applications. 
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Figure 1 | Schematic of surface-assisted bottom-up synthesis of 
atomically precise AGNRs from DBBA. a) side view of a ball- and- stick 
model of DBBA monomer, whose chemical structure is inserted; 
b) 7-AGNR; c) 14-AGNR; d) 21-AGNR. The number before "AGNR" 
indicates its width, that is, the number of carbon-pair lines in the direction 
perpendicular to AGNR axis, as sampled in Fig. Id. In Fig. lb, the yellow 
hexagon highlights an edged benzene-ring, the red arrows for the location 
between two benzene-rings edges. 

Results 

To provide a general characterization of the prepared AGNRs/ 
Ag(lll) sample, ex situ Raman measurements were carried out at 
ambient atmosphere. Fig. 2a shows an average Raman spectrum 
from ten spots. Six representative modes in the range from 
—350 cm" 1 to —1700 cm" 1 are observed and identified as G 
(—1610 cm" 1 ), D (—1350 cm" 1 ), the edge-phonon modes 
(—1220 cm" 1 and —1260 cm" 1 ) and radial-breathing-like-mode 
(RBLM) (-396 cm" 1 ) as well as its third overtone RBLM3 
(—960 cm" 1 ) 12 . Remarkably, the intensity ratio between modes D 
and G (Iv/Iq) is —0.90, which is about twice that of AGNRs on 
Au(lll) 7 . (also refer to Fig. SI) For polycrystalline graphite, the 
J D /J G ratio correlates with the effective crystallite size (La) as La 2 
for La < 2 nm: due to increasing fraction of sp 2 carbons, which 
contribute to the G band but not to the D band 1314 . According to 
Ref 14, the I D /I G ratio is -0.3 for 7-AGNR (0.74 nm) and -1.2 for 
14-AGNR. Thus, the larger I D /I G ratio on AGNRs/ Ag(l 11) suggests 
the coexistence of wider AGNRs such as 14-AGNRs or 21 -AGNRs 
with 7-AGNRs. Besides these six modes, the spectrum exhibits a 
sharp peak at —2000 cm" 1 , as highlighted in the middle inset, ori- 
ginating from the localized width-independent mode generated by 
the vibration of the atoms at the edge of AGNRs 15 . Notably, only 
when the AGNRs are at least partially H-free, that is, when C=C 
bonds exist at their armchair edges, can the localized mode be 
observed 15 . This gives evidence for the dehydrogenation at the head 
sites and the formation of AGNRs with multiple integer widths of 
14-, 21-AGNRs and so on. 

The STM image inserted in the upper- right corner in Fig. 2a 
clearly shows the randomly distributed AGNRs with varying widths 
on Ag(lll). The green rectangles highlight one 7-14-7-AGNR and 
one 14-21 -AGNR heterojunction, where the n- conjugation at the 
junction interface is preserved. The line profiles taken perpendicu- 
larly across 7-, 14-, 21-, and 28- AGNRs axes reveal the corres- 
ponding widths of —1.5 nm, —2.4 nm, —3.3 nm, and —4.2 nm 
respectively, and an apparent height of —0.19 nm, as shown in 
Fig. 2b. These values vary slightly depending on tip conditions. For 
such AGNRs the widths measured are larger than those from the 



corresponding atomic structure models. For example, the model of 
7-AGNR (Fig. lb) shows a smaller width of 0.741 nm. This could be 
attributed to the finite size of STM tip. The highly- resolved STM 
image taken at 77 K in Fig. 2c shows ESPs in 7- and 21-AGNRs, 
which arise from the intervalley electron backscattering at both arm- 
chair edges 10 ' 11 . 

Discussion 

We first discuss the two Raman peaks located at —2608 cm" 1 and 
—2690 cm" 1 in Fig. 2. A previous report shows Raman 2D peak 
splitting upon deposition of Au onto single layer graphene, which 
come from the graphene with and without Au, and is not the case in 
this work 16 . Recent theoretical and experimental investigations reveal 
that the Raman 2D peak of single layer graphene splits into two 
peaks under homogeneous uniaxial strain and both redshift signifi- 
cantly depending on the direction and magnitude of the strain 
applied 17 " 21 . In particular, the shift rates are up to —67.8 cm"V% 
and —26.0 cm" 7% for strain applied along the zigzag directions 17 . 
Accordingly, the average strain in the direction perpendicular to the 
AGNR axis is at least 2%. In this work, the mismatch between 
AGNRs and Ag(l 1 1) may induce strain. However, we do not observe 
obvious deformation of AGNRs along its zigzag direction. Further- 
more, uniaxial strain-induced redshift of the peak at —2690 cm" 1 is 
not observed. We also note that the 2D peak of GNRs formed by CNT 
unzipping is observed at —2659 cm" 1 6 . Thus, strain-induced Raman 
2D peak splitting is highly unlikely. We performed comparative 
Raman measurement on AGNRs of different widths. The Raman 
measurements on 7-AGNRs on Au(lll) is shown in Fig. SI, where 
one peak is observed at —2603 cm" 1 accompanied by a shoulder at 
—2685 cm" 1 . On Au surface, 7-AGNRs dominate, and the Raman 
peak at —2603 cm" 1 is attributed to 7-AGNRs. The shoulder at 
2685 cm" 1 appears to have the same origin as the peak at 
-2690 cm" 1 for AGNRs on Ag(lll). Wider AGNRs such as 14- 
and 21-AGNRs are observed on Ag(lll) compared to that on 
Au(lll). Hence, from our empirical observations, we suggest that 
the peak at —2690 cm" 1 is attributed to wider AGNRs. Further cal- 
culations are required to support this postulate. 

Fig. 3 displays STM images of ESPs in 7-, 14-, and 21-AGNRs (left 
column) as well as the corresponding DFT simulations (right col- 
umn), respectively. Figs. 3a and 3c are enlarged from Fig. 2c. ESPs in 
7-AGNR form localized ring-like shapes with a (^3 X V3)R30° unit 
cell, in good agreement with the calculated charge density plots of the 
conduction bands of 7-AGNR (filled red part in left panel in Fig. 3a) 7 
and with our DFT simulations (right panel in Fig. 3a). Such ESPs are 
attributed to the electronic states in the vicinity of the armchair edge 
where the incident and reflected waves interfere with each other due 
to intervalley scattering 1011 . A piece of ball and stick model of 7- 
AGNR is superimposed in the STM image of Fig. 3a, with the red 
dots highlighting each centre of the ESPs. This suggests that ESPs 
centres at both armchair edges of 7-AGNR are located not at but 
between the benzene-ring edges. One edged benzene-ring and the 
location between two adjacent benzene-ring edges are highlighted in 
Fig. lb by a yellow hexagon and red arrows, respectively. ESPs in 21- 
AGNR are similar to those in 7-AGNR, in a localized ring-like shape 
with a (^3 X y3)R30° unit cell, as marked by a red parallelogram. It is 
worth noting that there is a dislocation line parallel to the 21-AGNR 
axis, as highlighted by the green dash-dotted line. Given that 21- 
AGNR has the same lateral mirror symmetry as 7-AGNR and that 
the armchair edges have the same scattering effect on the incident 
electrons, we extend the above-mentioned calculated charge density 
plots of the conduction bands from 7-AGNR to 21-AGNR. The red 
(blue) dots indicate the ESP centres due to electron backscattering at 
the lower (upper) edge, as shown in the overlaid 21-AGNR ball and 
stick model in the left panel in Fig. 3c. Using geometrical considera- 
tions, ESPs from opposite armchair edges will be misaligned in the 
direction 60° off the lateral direction of 21-AGNR (parallel to the ^3 
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Figure 2 | a) Raman spectrum on AGNRs/Ag( 1 1 1 ) under ambient conditions. Besides mode D (1350 cm -1 ) and mode G (1610 cm" 1 ), some other 
vibration modes appear due to the finite width and low symmetry of the nanoribbons, such as RBLM (396 cm" 1 ), RBLM3 (961 cm" 1 ). The peak at 
2000 cm" 1 as highlighted in the middle inset arises from the vibration mode of the C-atoms at the edge of AGNRs. The inserted typical STM image (70 nm 
X 40 nm, V T = —2.0 V) at the upper-right corner shows straight AGNRs with varying widths, b) line-profiles taken perpendicularly across the 7-, 14-, 
21-, and 28-AGNRs, respectively, c) STM image (8 nm X 10 nm, V T = —0.5 V) showing the electron scattering patterns in 7- (lower) and 21 -AGNRs 
(upper). The 21-AGNR is imperfect at the lower end. 



direction as indicated by the red parallelogram), which will result in 
the observed dislocation. This is consistent with the DFT simulations 
in the right panel in Fig. 3c. The dislocation shifts to a lower position 
at the left end of the AGNR because of defects (see Fig. 2c, the lower 
end of 21-AGNR). 

ESPs in 14- AGNR are significantly different, as shown in the left 
panel in Fig. 3b. One-dimensional (ID) delocalized striped patterns 
were observed parallel to its axis with equal- spacing of —0.38 nm, in 
good agreement with the Fermi wavelength of graphene of 3a/2 = 
0.37 nm, where a = 0.246 nm is the lattice constant of graphene. 
Such ESPs arise from the intervalley electron backscattering at both 
armchair edges of metallic AGNRs 3 and depend not on the electron 
energy but only on the edge structure, as observed at the armchair 
edge of epitaxial graphene on SiC(0001) 22 . A piece of ball and stick 
model of 14-AGNR is overlaid to show the registration between the 
ESPs and C chains. This model matches perfectly with the predic- 
tion 3 : the probability density for a state near the Dirac singularity 



disappears at ]=3i (i.e. 3, 6, 9, 12 ), where i is an integer, and J is 

introduced to mark the position of carbon-pair lines in AGNRs, as 
shown in Fig. 3b. The good match with the DFT calculations in the 
right panel of Fig. 3b confirms the ID metallic nature of 14-AGNR 
due to quantum confinement effects 3 . 

In order to acquire more detailed information about the electronic 
structures of AGNRs, dl/dV measurements were performed at 4.7 K. 
These conductance measurements are directly proportional to the 
LDOS within the Tersoff-Harmann approximation 23 " 26 . The dl/dV 
spectra in Figs. 4, 5, and S2 reveal the electronic properties of such 
AGNRs. Fig. 4 summarizes the typical dl/dV spectra from the centres 
of three neighbouring AGNRs, which are shown in the inserted STM 
image. The bottom spectra taken from clean Ag(lll) terrace clearly 
show the surface states at —66 mV below Fermi energy level (E F , 0 V 
in the spectrum), as marked by a short cyan solid line, verifying the 
good tip conditions. The position of AGNR E F as indicated by the 
green dash-dotted line is shifted toward their unoccupied bands. This 
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Figure 3 | ESPs (left column, a: 3.5 nm X 2 nm, V T = -0.5 V; b: 5 nm X 
3 nm, V T = -0.5 V; c: 5.5 nm X 3.5 nm, V T = -0.5 V) and 
corresponding DFT simulations (right column) in 7-, 14-, and 21- 
AGNRs. Parameter J in panel b is an integer that marks the position of 
carbon-pair lines in AGNRs. The red parallelogram in Fig 3c marks the (^3 
X y3)R30° unit cell, the green dash-dotted line for dislocation, which shifts 
due to the defects at left end (see Fig. 2c). Red and blue dots: ESP centres. 

is attributed to the smaller work function of Ag(l 1 1) (—4.4 eV) and a 
hybridization between Ag d band and graphene 2p 27 - 28 . The spectra 

exhibit a sharp peak at 240 mV for 21-AGNR, a broad peak at 

400 mV for 14-AGNR and a broader peak at ~ -800 mV for 7- 

AGNR, as marked by green arrows, respectively. 

Fig. 5 shows dl/dV spectra acquired perpendicularly across 
AGNRs. In each panel, the black spectra were taken from the adja- 
cent Ag(l 1 1) surface, while the red (blue) ones are from the AGNRs 
edges (centre) at equal spacing. The inserted differential STM images 
show the typical AGNRs respectively. For 7-AGNR, dl/dV spectra 
exhibit a broad peak at 850 mV at the edges and is almost fea- 
tureless in the centre, as shown in Fig. 5a. Taking the substrate 
influence into account 28 , it is in consistence with the recent report 
by Ruffieux et al. 29 that the conduction band minimum of Au(lll)- 
supported 7-AGNR is located at 1.6 eV above E F at its armchair 
edges. It is also in good agreement with the recent report by Koch 
M. that the HOMO and LUMO of Au(l Unsupported 7-AGNRs 
delocalized along their armchair edges with an energy gap of 
2.7 eV 30 . Detailed inspection reveals that the unoccupied bands of 
7-AGNR is shifted below E F and thus becomes partially occupied, 
indicating a weak chemisorption nature 31 . On the contrary, while the 
featureless spectra are partially occupied at the edges, those in centre 

remain unoccupied and exhibit a broad (sharp) peak at 380 mV 

( 240 mV) for 14-AGNR (21-AGNR), as highlighted by solid 

green lines in Fig. 5 respectively. The spatially resolved spectra indi- 
cate that the electronic structures are heavily influenced when 
AGNRs adsorbed on Ag(lll). For 0 V < V T < 1.0 V, the LDOS 
of all the AGNRs are finite and similar to that of exposed Ag( 111) (for 
the case of 21-AGNR see Fig. S2). This suggests that in the energy gap 
of AGNRs the local density of states from underlying Ag(l 1 1) could 




Ag(111) 
_i . i . i i i 

-900 -600 -300 0 
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Figure 4 | STS data showing the unoccupied states of clean Ag(lll) and 
centers of 7-, 14- and 21-AGNRs, as marked on the inserted STM image 
(7 nm X 8 nm, V T = 1.5 V) as 1, 2 and 3 respectively, clearly show that 
electronic properties of AGNRs are significantly related to their widths. 

Offset vertically for clarity. The green dash- dotted line indicates E F and the 
short cyan solid line for the surface states of clean Ag( 1 1 1) at —66 mV. The 
green arrows indicate the typical peaks in STS spectra. 

not be excluded in STS measurements. Due to the charge transfer 
between AGNRs and Ag(l 1 1), the corresponding valence band max- 
imum of AGNRs shifts out of the measurement range and could not 
be detected with our experimental settings. Applying Ruffieux's 
method 29 , the energy gap of semiconducting 7-AGNRs is deduced 
to be over 1.8 eV, in good agreement with recent experimental 
reports 29 ' 30,32 ' 33 . For 21-AGNRs, the energy gap is over 1.2 eV. For 
metallic 14- AGNRs, the observed peak at —380 meV might be the 
first Van Hove singularity (VHS) at the unoccupied states, which is 
broadened due to the hybridization with Ag(lll) 28 ' 34 . The first VHS 
gap of metallic 14-AGNR is calculated to be close to 2 eV 35 , and is 
also out of the measurement range. Further experiments and simula- 
tions are required to accurately measure the energy gaps. 

For heterojunctions such as 7-14-AGNRs, 14-21-AGNRs, the 
spectra taken along the ribbons exhibit the previously mentioned 
characteristics of each AGNR when measured in their respective 
regions (Fig. 4). The change in spectra across the boundary (< 
1 nm) is very sharp, having no transitional spectra. This is similar 
to the case of metallic intramolecular heterojunctions in CNT 36 , 
confirming that the electronic properties of AGNRs are significantly 
width-dependent. 

Why do AGNRs with varying widths form on Ag( 1 1 1 ) ? To answer 
this question, a submonolayer DBBA on Ag(lll) annealed at 
—450 K for 30 min was investigated. Unexpectedly, instead of 
anthracene polymers, isolated planar products are observed, as 
shown in Fig. 6, indicating that dehydrogenation has taken place 
for DBBA monomers on Ag(lll) at this stage 37 . The carbon back- 
bone of an isolated molecular product is inserted at the upper-right 
corner of Fig. 6a. More theoretical simulations are required to deter- 
mine how many hydrogen atoms dissociate from one DBBA mole- 
cule and its final configuration. Two such structural units are 
tentatively superimposed on the high resolution STM image in 
Fig. 6d. Fig. 6b shows the line profile taken along the white line in 
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Figure 5 | STS spectra along the lines across 7-, 14- and 21-AGNRs perpendicularly, respectively. Offset vertically for clarity. The green dash-dotted 
lines indicate the Fermi level, the green lines for the main peaks. 



Fig. 6a, and Fig. 6c shows the corresponding fast Fourier transform 
(FFT) pattern from Fig. 6a, revealing that isolated planar molecular 
products form a hexagonal arrangement with a nearest neighbour 
distance of 1.5 ±0.1 nm. This is much larger than that of — 1.15 nm 
for monolayer coronene on Ag(l 1 1) 38 , indicating a strong adsorbate- 
substrate interaction at 77 K. Some short rod-like features of two or 
three connected intermediate products are observed in the area 
enclosed by white dash- dotted curve, indicating that such planar 
products act as intermediates for AGNRs. The growth mechanism 
of AGNRs on Ag(lll) is proposed in Fig. S3, which is significantly 
different with that on Au(lll) 7 . It may be related to the different 
catalytic reactivities of Au(l 1 1) and Ag(l 1 1), which induce different 
self-assembly behaviour of DBBA at room temperature, as shown in 
Fig. S4. 

In summary, we have succeeded in fabricating ultra-narrow atom- 
ically precise 7-, 14- and 21-AGNRs and observed their associated 
ESPs using LT-STM. The observed ESPs are in good agreement with 
theoretical predictions. The dl/dV spectra reveal the width-dependent 



electronic properties of the AGNRs. For the next step of our invest- 
igation, we plan to investigate the electronic properties of such 
AGNRs on insulating thin films. We expect that these initial investi- 
gations on the electronic properties of such AGNRs give us a better 
understanding of this extraordinary class of materials and help ad- 
vance theoretical predictions on systems such as intraribbon 
quantum dots 39 , GNR switches 40 and magnetic devices based on spe- 
cific GNR edge states 41 . 

Methods 

Synthesis of AGNRs. The clean Ag(lll) surface was achieved in situ by repeated 
cycles of Ar+ sputtering and subsequent annealing to —750 K. DBBA (97%, St- Jean 
Photochemical Inc.) was evaporated from a K-cell (MBE-Komponenten, Germany) 
at a deposition rate of —0.01 ML/min (ML: monolayer) at normal incidence onto 
Ag(l 1 1) at room temperature. After deposition, the sample was annealed at —450 K 
and —650 K for for 30 mins, respectively. 

STM/STS measurements. Sample preparation and STM measurements were carried 
out in a custom-built multi-chamber ultra-high vacuum (UHV) system housing an 




Figure 6 | Dehydrogenation of DBBA on Ag(lll) at —450 K. a) STM image (V T = — 1.0 V) showing that most of the products are isolated in a 
hexagonal arrangement, as revealed by the corresponding FFT pattern in panel c; while some are connected together as short rods, as highlighted by the 
white dash-dotted curve. The carbon backbone of the isolated product is inserted at the upper right corner, b) the line-profile taken along the white line in 
panel a, showing an average intermolecular distance of 1.5 ±0.1 nm. d) zoomed-in STM image (V T = —0.01 V) showing the isolated product lies flat on 
Ag(lll). Two molecular unit structures are superimposed. Such planar products are believed to be the intermediates for AGNRs in varying widths. 
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Omicron LT-STM in the analysis chamber with a base pressure better than 

6.0 X 10" 11 mbar. All the STM images were recorded in constant current mode using 

chemically etched tungsten (W) tips. To reduce the thermal drift and noise level, 

differential conductance (dl/dV) measurements were performed at —4.7 K in an 

open feedback loop mode using the lock- in technique with frequency 625 Hz and ac 

modulation of 20 meV and with set-point parameters of V T = 100 mV and 

I = 200 pA. 

Raman measurements. Raman spectra were recorded with a Renishaw 1000 Raman 
Microscope at ambient. The excitation source was a 514 nm laser with power below 
20 mW to avoid laser induced surface heating. The laser spot size was around 10 um. 

Simulation of the ESPs. Simulations of electron scattering patterns (ESPs) of pristine 
graphene nanoribbon were obtained by performing ab-initio pseudopotential density 
functional theory calculations. In Tersoff and Hamann approximation 17 , the 
tunneling current density between the STM tip and the sample is proportional to the 
local density of states (LDOS) 

JEp-Vbias 

<&i, R is the wave function of pristine graphene nanoribbon. When a small bias voltage 
Vbias is applied, the numerical integration is carried out on the associated eigenvalues 
within this energy range below the Fermi level. This provides occupied (unoccupied) 
electronic states when a negative (positive) bias voltage is applied. The exchange- 
correlation functionals are described by the generalized gradient approximation 
(GGA) derived Perdew-Burke-Ernzerhof (PBE) approach 18 . The calculations were 
performed by using the SIESTA code 19 , together with a double-^ polarized basis-set 20 , 
a real space mesh cutoff of 400 Ryd and norm conversing pseudopotential. The slab 
models include a vacuum thickness of 10 A to reduce the influence of periodic images 
along the direction perpendicular to the GNR plane. The total energies converged well 
with respect to the sufficiently large k-point sampling of 90 X 1 X 1 in the Brillouin 
zone. The electronic structure calculations were carried out using fully relaxed 
structures where the structural optimization on the GNR atomic degrees of freedom 
was achieved until the residual forces on each atom are smaller than 0.01 eV/A. 
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